P.B.5818 - Palenilaan 2 
2260 HV Rijswfjk (ZH) 
2? +31 70 340 2040 
TX 31651 epo nl 
FAX +31 70 340 3016 



EuropSisches 
Patentamt 



Zwelgsteie 
in Oen Haag 
Recherchen- 
abteflunq 



European 
Patent Office 

Branch al 
The Hague 
Search 
^division 



Offlc uropeen 
ctesbrev ts 



Departement a 
La Haye 
Division de la 
fftchfifchfi 



r Holt, Michael 
Texas Instruments Ltd., 
EPD MS/13, 
800 Pavilion Drive 
Northampton Business Park, 
Northampton NN4 7YL 
GRANDE BRETAGNE 



DalunVDale 

22.04.03 



L 


j 1 — 1 


ZeichervTtofyfttil. 

TI-33782 EP 


Anmeidung NriApplication NoVDemande n* JPalenI Nr ./Patent NoVBrevel n*. 

02293152.1-1237- 


AnmeWer/Afmlicant/Dernandeur/Palerri 

Texas Instruments Incorporat 


:ed, et al 


COMMUNICATION 



The European Patent Office herewith transmits as an enclosure the European search report for the 
above-mentioned European patent application. 

If applicable, copies of the documents cited in the European search report are attached. 

Q Additional set(s) of copies of the documents cited in the European search report is (are) enclosed 
as well. 

The following specifications given by the applicant have been approved by the Search Division: 
□ abstract □ title 

DQ The abstract was modified by the Search Division and the definitive text is attached to this 
communication. 

The following figure will be published together with the abstract 4 



REFUND OF THE SEARCH FEE 



DOCKETED 

29 APR 2003 

IMPETUS 



Received 

23 :::-3 
Impetus 



If applicable under Article 10 Rules relating to fees, a separate communication 
from the Receiving Section on the refund of the search fee will be sent later. 




European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Numbar 

EP 02 29 3152 



DOCUMENTS CONSIDERED TO BE RELEVANT 




Category 


Citation of document with indication, where appropriate, 
of relevant passages 


Relevant 
to claim 


CLASSIFICATION OF THE 
APPLICATION (lntCI.7) 


X 


WO 02 39689 A (ERICSSON) 
16 May 2002 (2002-05-16) 

* page 3, line 24 - page 4, line 2 * 

* page 4, line 10 - line 14 * 

* page 4, line 25 - line 27 * 

* page 12, line 7 - line 9 * 

* page 15, line 1 - line 4 * 

* page 15, line 10 - page 16, line 4 * 


1-10 


H04L25/03 


X 


W0 01 60002 A (ERICSSON) 

16 August 2001 (2001-08-16) 

* page 8, line 27 - page 9, line 3 * 


1-10 




X 


W0 01 03308 A (FRAUNHOFER) 
11 January 2001 (2001-01-11) 

* page 1 , paragraph 2 * 

* page 2, paragraph 3 * 


1,6 




X 


W0 01 78239 A (QUALC0MM) 
18 October 2001 (2001-10-18) 


1,6 






* page 4, line 35 - page 5, line 9 * 




TECHNICAL FIELDS 
SEARCHED (lntCl.7) 


X 


FETTWEIS ET AL. : "Reduced-complexity 
V1terb1 detector architectures for partial 
response signalling" 

IEEE GLOBAL TELECOMMUNICATIONS CONFERENCE, 

14 - 16 November 1995, pages 559-563, 
XP000621547 
New York, US 
* figures 2-4 * 


1,6 


H04L 
H03M 


X 


EP 1 058 392 A (MOTOROLA) 
6 December 2000 (2000-12-06) 
* page 2, line 23 - line 31 * 

-/- 


1,6 




The present search report has been drawn up for afl claims 







Raoo of search 

THE HAGUE 



Data of comp^tion at tha March 

3 April 2003 



Examiner 

Scrlven, P 



CATEGORY OF CITED DOCUMENTS 



X : particularly relevant it taken alone 

Y : particularly relevant if combined with another 

document of the same category 
A : technological background 
O : non-written disclosure 
P : Intermediate document 



T : theory or principle underlying the Invention 
E : earlier patent document, but published on. or 

after the Tiling date 
D : document died In the application 
L : document cited for other reasons 

& : member of the same patent ramify, corresponding 
document 



page 1 of 2 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application Number 

EP 02 29 3152 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation ol document with Indication, where appropriate. 
of relevant passages 



PAGE, CHAU: "Improved architectures for 
the add-compare-select operation 1n long 
contstralnt length Viterbl decoding" 
IEEE JOURNAL OF SOLID-STATE CIRCUITS, 
vol. 33, no. 1, January 1998 (1998-01), 
pages 151-155, XP000766629 
Mew York, US 
ISSN: 0018-9200 

* page 151, left-hand column, paragraph 3 
- right-hand column, paragraph 1 * 

EP 0 786 872 A (LUCENT) 
30 July 1997 (1997-07-30) 

* page 5, line 31 - line 55 * 

SHEEN, STUBER: "Error probability of 
reduced-state sequence estimation for 
trellis-coded modulation on Intersymbol 
Interference channels" 
IEEE GLOBAL TELECOMMUNICATIONS CONFERENCE, 

2-5 December 1991, pages 1159-1164, 
XP010042585 
New York, US 

* figure 3 * 

PR0AKIS J G: "Digital Communications" 
DIGITAL COMMUNICATIONS; MCGRAW-HILL, 
1995, pages 
249-254 , 468-470 , 621-627 , 649-654 , 69-8-702 , E 
TC-, XP002202678 
New York, US 
ISBN: 0-07-113814-5 

* page 249 - page 254 * 



The present search report has been drawn up for all claims 



Place ol soarch 

THE HAGUE 



Relevant 
to claim 



1,6 



1-10 



2-5,7-10 



1,6 



Data of completion of the search 

3 April 2003 



CLASStFlCATION OF THE 
APPLICATION <lnLCK7) 



TECHNICAL FIELDS 
SEARCHED flntCI.7) 



Examiner 

Scrlven, P 



CATEGORY OF CITED DOCUMENTS 



X : particularly relevant if taken alone 

Y : particularly relevant if combined with another 

document of the same category 
A : teennotogteai background 
O : non-written disclosure 
P : Intermediate document 



T : theory or principle underlying the Invention 
E : earlier patent document, but published on, or 

after the Wing date 
D : document cited In the appfication 
L : document cited tor other reasons 

A : member of the same patent family, corresponding 
document 



page 2 of 2 



ANNEX TO THE EUROPEAN SEARCH REPORT 
ON EUROPEAN PATENT APPLICATION NO. 



EP 02 29 3152 



This annex lists the patent family members relating to the patent documents cited in the above-mentioned European search report. 
The members are as contained in the European Patent Office EDP file on 

The European Patent Office is in no way liable for these particulars which are merely given tor the purpose of information. 



03-04-2003 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



W0 0239689 



16-05-2002 



AU 
WO 



1451102 A 
0239689 A2 



21-05-2002 
16-05-2002 



WO 0160002 


A 


16-08-2001 


AU 


3282601 A 


20-08-2001 








EP 


1256217 Al 


13-11-2002 








WO 


0160002 Al 


16-08-2001 


WO 0103308 


A 


11-01-2001 


WO 


0103308 Al 


11-01-2001 








AU 


5032299 A 


22-01-2001 








EP 


1192719 Al 


03-04-2002 


W0 0178239 


A 


18-10-2001 


AU 


5142801 A 


23-10-2001 








EP 


1269636 A2 


02-01-2003 








WO 


0178239 A2 


18-10-2001 


EP 1058392 


A 


06-12-2000 


EP 


1058392 Al 


06-12-2000 


EP 0786872 


A 


30-07-1997 


EP 


0786872 A2 


30-07-i997 








JP 


10004362 A 


06-01-1998 



Si For more details about this annex : see Official Journal of the European Patent Office, No. 12/82 



ABSTRACT / ZUSAMMENFASSUNG / ABREGE 



02293152.1 



A Reduced State Sequence Equalizer (RSSE) is implemented using a butterfly hardware 
accelerator (58) to increase the efficiency of computing branch metrics (24) and the ACS 
(add.compare and select) function (26). 



1 



(19) 



J 



(12) 



(43) Date of publication: 

30.07.1997 Bulletin 1997/31 

(21) Application number: 96309397.6 

(22) Date of filing: 23.12.1996 



Europdlsches Patentamt 
European Patent Office 
Office europ^en des brevets (11) EP 0 786 872 A2 

EUROPEAN PATENT APPLICATION 

(51) Int. CI. 6 : H03M 13/00 



(84) Designated Contracting States: 


(72) Inventor: Cesarl R.A. 


DE ESFRGBIT 


Murray Hill, New Jersey 07974-0636 (US) 


(30) Priority: 29.12.1995 US 9537 


(74) Representative: Johnston, Kenneth Graham et at 




Lucent Technologies (UK) Ltd, 


(71) Applicant: Lucent Technologies Inc. 


5 Mornington Road 


Murray Hill, New Jersey 07974-0636 (US) 


Woodford Green Essex, IG8 OTU (GB) 



(54) Viterbl decoder with reduced metric computation 
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computing branch metrics that are combinations of 
autocorrelation and cross correlation terms are con- 
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the autocorrelation and cross correlation terms at a 
symbol instant. Once calculated and stored, predeter- 
mined ones of the autocorrelation components and the 
cross correlation components are selected. The 
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tion components, or their inverse as predetermined, are 
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mined combinations of the stored components of auto- 
correlation and cross correction terms, or their inverse, 
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same symbol instant. All branch metrics associated with 
the symbol instant can be calculated in this manner. The 
components of the autocorrelation terms and cross cor- 
relation terms associated with the next symbol instant 
are calculated and stored, and the process of generat- 
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Technical Field 

5 This invention relates general!/ to viterbi decoders such as would be used in digital receivers, and particularly to a 
viterbi decoder with reduced metric computations. 

Background pf the Invention 

w Although telecommunications technology originally developed using copper wires as the transmission media, today 
there is much interest and commercial activity in wireless technology using free space as the transmission media. Wire- 
less technology frees the user from the constraint of needing a physical connection to the transmission media. 

Initial commercial activity developing wireless technology used analog transmission technology. In the interest of 
increasing efficiency of spectrum utilization, digital transmission systems have been developed. One digital standard 

is (Group Special Mobile, or GSM) has been adopted for use in Europe. A different standard (IS54) will likely be adopted 
for use in North America: either a time-division multiple access (TDM A) standard, or a competing code-division multiple 
access (CDMA) standard. Other parts of the world such as Japan are considering yet another digital standard, such as 
the Japan Digital Standard. 

Digital systems transmit information either as zeros or ones. While accurate transmission of information is generally 

20 not at issue, accurate detection of information received depends on being able to detect transmitted ones as ones upon 
being received, and being able to detect transmitted zeros as zero upon being received. Accurate detection is not a sim- 
. pie matter due to many factors. These factors include low received signal power at a receiver, movement of the trans- 
mitter relative to the receiver which introduces changing channel characteristics. Multi-path interference, or inter-symbol 
interference, occurs when reflections from various objects, including atmospheric layers, cause the presence of multiple 

25 transmission paths from the transmitter to the receiver, with varying transit times for a transmitted signal, causing sev- 
eral images of a transmitted signal to arrive at the receiver at different instants in time. 

Error correction techniques have been developed to increase accuracy of digital communication systems. A viterbi 
decoder is a maximum likelihood decoder that provides forward error correction. Viterbi decoders are used to decode 
a sequence of encoded symbols, such as a bit stream. The bit stream can represent encoded information in a telecom- 

30 munications system. Such encoded information can be transmitted through various media with each bit (or set of bits) 
representing a symbol instant. In the decoding process, a viterbi decoder works back through a sequence of possible 
received bit sequences at each symbol instant to determine which one bit sequence was most likely to have been trans- 
mitted. The possible transition from a bit at one symbol instant, or state, to a bit at a next or a subsequent, symbol 
instant or state, is limited. Each possible transition from one state to a next state can he shown graphically and is 

35 defined as a branch. All possible branches constitute a trellis. A sequence of interconnected branches is defined as a 
path. Each state can transition only to a limited number (for example two or four) of next states upon receipt of the next 
bit (or set or bits) in the bit branch. Thus, some paths survive and other paths do not survive during the decoding proc- 
ess. By eliminating those transitions that are not permissible, computation efficiency can be achieved in determining the 
most likely path to survive. The viterbi decoder typically defines and calculates a branch metric associated with each 

40 branch and employs this branch metric to determine which paths survive and which paths do not survive. 

A branch metric is calculated at each symbol instant for each possible branch. Each path has an associated metric, 
accumulated cost, that is updated at each symbol instant. For each possible transmission, the accumulated cost for the 
next state is calculated as the extremum of the sum of the branch metric for the possible transmission and the path 
accumulated cost at the respective previous, or originating, state. Calculation at each symbol instant of the branch met- 

45 ric involves an add, compare and select operation for each state. This can result in computation-intensive operation. 

A need remains in the art for all efficient way to calculate the branch metrics that substantially reduces the number 
of calculations and concomitantly reduces the power and time consumed in the calculation of branch metrics. 

Summary pf the Invention 

50 

In accordance with the invention, a method and apparatus for generating branch metrics that are combinations of 
autocorrelation and cross correlation terms are constructed by first calculating and storing components of the autocor- 
relation and cross corr lation terms at a symbol instant. Once calculated and stored; predetermined ones of the auto- 
correlation components and th cross corr lation components are selected. The selected autocorrelation components 
55 and cross correlation components, or their inverse as predetermined, are combined to produce a branch metric. Other 
predetermined combinations of the stored components of autocorrelation and cross correlation terms, or their inverse, 
are combined to produce other branch metrics at the same symbol instant. All branch metrics associated with the sym- 
bol instant can be calculated in this manner. The components of the autocorrelation terms and cross correlation terms 
associated with the next symbol instant are calculated and stored, and the process of generating the branch metrics 
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repeated. 

Brief Description of the Drawings 

5 The invention wiil now be described by way of example in which : 

Figure 1 is a block diagram of a transceiver, in accordance with an illustrative embodiment of the present invention; 
Figure 2 is a diagram showing three symbol instants and the possible branches among the individual states; 
Figure 3 is a graphical representation of an odd symbol phase constellation for n/4- shifted differentially encoded 
10 quadrature phase shift keying; 

Figure 4 is a graphical representation of an even symbol phase constellation for n/4- shifted differentially encoded 
quadrature phase shift keying; and 

Figure 5 is a block diagram of a circuit for calculating branch metrics in a viterbi decoder, using reduced computa- 
tions in accordance with the present invention; and 
is Figure 6 is a graphical comparison of the calculation of branch metrics in accordance with the present invention, as 
compared to the prior art. 

Detailed Descript ion 

20 Transceiver 1 0 shown in the block diagram of Figure 1 is an illustrative embodiment of the invention. Transceiver 1 0 
is a communication system or portion of a communication system such as a mobile digital cellular telephone. Trans- 
ceiver 10 is comprised of transmitter 12 and receiver 14. 

As shown in Figure 1 , a microphone 121 for converting a voice signal into an electrical signal is coupled to a voco- 
der 122. The vocoder is coupled to the channel encoder 123. which includes reordering and partitioning 1231, inter- 
ns leaving and formatting 1 232, and symbol converter 1 234. The output of the channel encoder is coupled to a transmitter 
filter 124. The filter output is coupled to a continuous balanced modulator 125, which in turn is coupled to another filter 
126. Filter 126 is coupled to an up converter and power amplifier 127 which is turn is coupled to an antenna 128. 

On the receiver side, antenna 141 is coupled to a low noise amplifier and down conversion 142, which is coupled 
to f Htering 143. Filtering 143 is coupled to automatic gain control 144 which conditions the signal for input to the Viterbi 
so decoder 145. The Viterbi decoder is coupled to a channel decoder 146 which includes deinterleaving, deformatting, 
deordering and departitioning 1461. The channel decoder is coupled to a vocoder 147. The vocoder is coupled to a 
speaker 148. 

Figure 2 shows a trellis diagram useful in understanding the branch metrics used in a Viterbi decoder. A Viterbi 
decoder reconstructs a transmitted bit stream from a received bit stream. Perfect digital signal symbols in a sequence 

35 are transmitted. However, the received signal is distorted due to noise and interference. As the mobile digital cellular 
telephone moves through its environment, the characteristics of the channel over which the signals are transmitted and 
received change. At each state or symbol instant (shown horizontally), there are a number of states (shown vertically) 
that range from zero up to 4 C " , -1 where C is the constraint length. These 4 C I states are referred to as individual states. 
The portion of a trellis diagram shown in Figure 2 is constructed one symbol instant at a time, as is known in the art. 

40 The trellis entries at each symbol instant comprise a cell. At each symbol instant, each individual state of the illustrative 
embodiment shown in Figure 2 can only transition to four possible individual states at the next symbol instant. Concom- 
itantly, each individual state at a symbol instant has only four possible previous individual states from which a transition 
to it can originate. Other decoders could have more or less than four transitions from a previous state and to a next 
state. 

45 A branch metric is calculated at each symbol instant for each possible transition from an individual state at one 
symbol instant, known as a source state, to an individual state at the next symbol instant. Various methods for calculat- 
ing branch metrics are known in art. The branch metric of all branches transitioning to a given next individual state are 
calculated, then added to the accumulated cost of the respective originating individual state. This results in four poten- 
tial accumulated cost sums. A comparison is made of the four potential accumulated cost sums. An extremum branch 

so metric, such as the least of the four sums, is selected as the next state accumulated cost A transition corresponding to 
the least sum is the most likely transition from the four possible originating individual states to the given next individual 
state. The surviving branch data of the most likely transition is stored forming the trellis of Figure 2 as the surviving 
branch to the given next individual state as is known in the art. The extremum branch metric is the accumulated cost of 
the next individual state. This process is repeated for each of th next individual states and repeated also for each sym- 

55 bol instant in the sequence. 

Accumulated cost and the source state for th selected branch can be stored in an array. When the trellis is com- 
plete, a decoded symbol at each symbol instant can be obtained by a traceback operation. Starting with the individual 
state in the last symbol instant having the minimum associated accumulated cost, a single trace back through the trellis 
will provide a decoded symbol at each symbol instant. At each symbol instant of this traceback operation, the source 
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15 



state is translated into a decoded symbol for that symbol instant, as is known in the art. 

A viterbi decoder is used with a partial log likelihood estimator to determine the most likely symbol sequence taken 
from a set of possible sequences over an interval of time. Minimizing the integral squared error between the received 
sequence and the expected sequence results in the most likely sequence to have been transmitted through a channel. 
The predicted squared error or log-likelihood between the received signal and an expected sequence (or complex func- 
tions is the auto correlation of the predicted error functions, multiplied by a constant. The constant may be ignored when 
evaluating relative likelihoods. 



log likclihood m = 



N-i l-I 

X y(n) - X s ^)m . Mn-ri) 



N-1 l-I 
n=0 i-0 



(I) 



where: 



20 



25 



30 



y(n) is a received sample sequence 

S(i) is the predicted symbol at symbol instant or symbol index i 

h is a channel estimate, estimated at the symbol rate (with index i) 

iis the symbol index 

I is the number of symbols in the interval 

n is l/r 

N is the number of samples in a sequence, N=l/r 
m is the metric, or expected value index 
r is the number of samples per symbol 
* represents a complex conjugate. 

Recognizing and taking advantage of conjugate symmetry, and substituting the index j for the index i in the conjugate 
error factor, the log-likelihood can be rewritten as 



35 



40 



log likclihoud m = 
«r- 1 < N-l l-l 

= X yU) 2 -2RcXyOO < XS(i) 



. h(ft-rif 
kji-in-i r n 

+2R*x X I S(i)"S(j)1 * h < n - ri > h < n ~ r i> 



(2) 



45 



50 



55 



The expected value (the right term in the right bracket of equation (1)) can be considered a sum of partial responses, 
each of which is equal to S(i) m multiplied by the channel response h(n) and placed in a time location r * i. For a channel 
filter, L symbols long, where L is the number of symbols in the channel impulse response, L partial responses overlap. 
The received signal is processed r samples at a time. The total squared error, which is proportional to the log-likelihood, 
is the sum of the partial windowed squared error values. For a 4-ary case, each symbol can take on one of four possible 
values. With three symbols affecting each one-symbol window, each window can have 64 possible partial log-likeli- 
hoods. For a rapidly changing channel, it may be necessary to update the channel estimate, h(n), once per symbol. 

Equation (2) represents the integral of the signal energy minus the difference between the cross correlation of the 
expected sequence. The three terms of equation (2) represent energy, cross-correlation, and autocorrelation respec- 
tively. Terms not dependent on S(i) m need not be included when calculating the relative probabilities of the symbol 
sequences. Such terms include the energy term and the inner products of the autocorrelation term where i=j. It is only 
necessary to find the sequence which maximizes the absolute value of the difference of the last two terms in the above 
equation. 

The equation for evaluating the partial log-likelihood for a three-tap channel filter, 16-state case are: 
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log likelihood (m) = S(i-I) m . y(n) * . h ^(k+l) + S(i) m . y(n)* . h n (k) 
+S(i+l) m . y(n)* . h^Oc-Q + [S(i-I) . S(i)] m . h^k+l) . h n (k+l) + [Stf-I)' . S(i+I)j m . h^W . h (n+0 (k) 

+ [SOr.S(W)] m .h n (k-l)*.h n+l (k-l) 

5 

where: 

•Vq W is tne P as t channel estimate, 
h( 0) (K) is the present channel estimate, and 
10 n( n+ ij (k) is an approximation of a future channel estimate. 

When performing the initial timing and channel estimate, a known symbol sequence is compared with the received sig- 
nal to obtain the most likely channel impulse response. For synchronous detection with the differential quadrature 
phase shift keying (DQPSK) channel, the known preamble sequence of Differential symbols should be first converted 

is into a set of absolute phase values from which an expected value is generated for use in the channel estimation proc- 
ess. The first symbol in the preamble used for the channel estimate should be an even number of symbol back from the 
first unknown data symbol, and the value of this first phase value prior to the first known symbol used for the channel 
estimate should be arbitrarily set to zero phase. The DQPSK phase differential symbol values following the first symbol 
preamble used in the channel estimation should then be converted to absolute phases accordingly. The phase offset of 

20 the first received symbol can then be expected to be one of the even indexed phases. 

[2n n/8; n even] 

on the unit circle (when using the channel estimate derived from the above-mentioned preample sequence). The 
25 expected complex symbol value for the first unknown symbol is then one of the four complex values; 1+Oj 0-1 j. -1+0j. 
and 0+1 j. 

Odd indexed symbols, are expected to be one of four odd phases on the unit circle 

[2k n/8; n odd] 

30 

By rotating the odd indexed expect partial responses by -rt/4, the MLSE operation remains identical for both even and 
odd indexed symbols, and eliminates the need for performing a scalar operation for each state in a symbol interval. The 
- ti/4 phase rotation requires calculating the sum and the difference of the real and imaginary value of the expected val- 
ues and multiplying then by V2, and re-ordering these scaled sums and differences in four possible arrangements with 
35 the appropriate sign changes. Table I summarizes the arithmetic for eight phase rotations on the unit circle. 



TABLE I 



0 in radians 


ReteHRcA+jlmA)] 


Imfe^ReA + jlmA)] 


0 


Re A 


ImA 


n/4 


V2 (Rc A-lm A) 


V2 (Re A+lm A) 


2rt/4 


-ImA 


Re A 


3n/4 


V2(-Re A-lm A) 


V2 (Re A+lm A) 


4n/4 


-Re A 


-ImA 


5n/4 


V2 (-Re A+lm A) 


V2(-Re A-lm A) 


6rt/4 


ImA 


-Re A 


7n/4 


V2 (Re A+lm A) 


V2 (-Rc A+lm A) 



The combination of appropriately chosen references greatly simplifies the calculation of the partial log-likelihoods. For 
55 example, th r samp! sum-of-products for each term is calculated only once. Each t rm, then reduced to a single com- 
plex value, assumes the real value or one of four rotations. This operation simply involves selecting either the real of 
imaginary valu of the complex number and multiplying it by + 1 or - 1 . The partial log-likelihoods are assigned as met- 
rics to a 16 state trellis shown in Figure 2. Four paths with unique metrics enter each state. Four paths with uniqu met- 
rics also originate at each state. For example, a channel having a multipath spread of three symbols values affect each 
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partial log-likelihood window. Viterbi decoding is applied to the trellis to determine the path with the largest accumulated 
partial log likelihood. Starting with the state in the last symbol instant of the interval having the largest accumulated par- 
tial log-likelihood, a path is traced back to reveal the most likely symbol sequence for the interval by decoding a symbol 
at each symbol instant in the traceback operation. The symbol sequence is then converted to differential symbol values 
for final DQPSK bit detection. 

A circuit useful in generating the 64 partial log likelihood values, which are referred to as incremental metrics of 
branch metrics, for each symbol instant employing equation (2) in accordance with the invention is shown in Figure 5. 
For a channel impulse response having three terms, that is a channel filter having three taps, and one sample per sym- 
bol, in equation (2) r=l, i=N and UN. 

The past and present channel estimates, as well as an approximation of a future channel estimate are provided at 
500, 502 and 504. At each symbol instant, the reconstructed symbol is created for a single integration window. The 64 
incremental metrics are developed from the integration window. The three taps of the channel impulse response are 
provided in each of blocks 500, 502 and 504. Both the in-phase I, and quadrature, Q. components of a different one of 
the channel impulse response taps provides an input to each rotate block. In a preferred embodiment, the odd indexed 
channel estimates are rotated back by k/4, in respective rotate blocks 506, 508 and 510, to provide an expected value 
at 0, 2rc/4, etc. Rotation is accomplished by substitution using TABLE I. Since only the odd channel estimates are 
rotated, at one symbol instant ROTATE blocks 506 and 510 will rotate channel estimates, and at previous and subse- 
quent symbol instants only ROTATE block 508 will rotate the channel estimate. 

Each rotate block generates in-phase and quadrature outputs of the partial response. Outputs 512 and 514 from 
rotate block 506 are provided to both complex conjugate generator 524 and complex multiply 526. Complex conjugate 
generator 524 produces the complex conjugate of its inputs 512 and 514 as its outputs 526 and 528 and provides the 
outputs as two inputs to complex multiply 530. The other two inputs to complex multiply 530 are outputs 516 and 518 
from rotate block 508. Complex multiply 530 computes the complex product of its inputs and provides in-phase and 
quadrature outputs 532 and 534, which are stored in registers 564 and 566. 

Outputs 516 and 518 from rotate block 508 are provided to both complex conjugate generator 536 and to complex 
multiply 538. The complex conjugate of inputs 51 6 and 51 8 are provided as outputs 540 and 542 from complex conju- 
gate generator 536 which are two inputs to complex multiply 544. The other two inputs to complex multiply 544 are out- 
puts 520 and 522 from rotate block 510. Complex multiply 544 computes the complex product of its inputs and provides 
in-phase and quadrature outputs 546 and 548 which are stored in registers 568 and 570. 

Outputs 520 and 522 from rotate block 510 are provided to both complex conjugate generator 550 and to complex 
multiply 552. The complex conjugate of inputs 520 and 522 are provided as outputs 554 and 556 from complex conju- 
gate generator 550 which are two inputs to complex multiply 558. The other two inputs to complex multiply 558 are out- 
puts 512 and 514 from rotate block 506. Complex multiply 558 computes the complex product of its inputs and provides 
in-phase and quadrature outputs 560 and 562 which are stored in registers 572 and 574. 

Complex multiplys 526, 538 and 552 each multiply the inputs from respective rotate blocks by sampled received 
signal y n \ and respectively store the complex product in registers 576 and 578, 580 and 582, and 584 and 586. 

The complex multiply function represented by complex multiply blocks 526. 530, 538, 544, 552 and 558 can be per- 
formed by a single multiplier. Similarly, other functions represented in Figure 5 need not be performed by separate hard- 
ware. 

State table control 588, shown in Table II, controls multiplexers 590-600 by select inputs 602-612. Each multiplexer 
output 61 4-624 is provided to a respective selector 626-636. Selectors 626-636 are also controlled by state table control 
588 by select input 638-648. Each selector 626-636 provides an output 650-660 to adder 662 which adds each of the 
inputs to produce output 664. 
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TABLE 1! 



The control of multiplexers 590-600 and selectors 626-636 is represented in Table 1 1. The first column is the metric 
25 counter. The second, third and fourth columns are the symbol indices. The six columns at the right represent the twelve 
controls 602-612 and 638-648. The I or Q in the first of the six right-most columns represents whether the in-phase or 
quadrature output from registers 576 or 578 is selected as the input to multiplexer 590. The sign preceding the entry in 
the first of the six right-most columns represents whether the selector 626 selects the output of the multiplexer taken as 
positive for a +, or inverted for a - . The other five of the six right-most columns each control one of the other muftplexer- 
30 selector combinations. 

The control information in the six right-most columns of Table II to control multiplexers 590-600 and selectors 626- 
636 can be stored in an array with each branch metric requiring only twelve control bits, two for each term of equation 
(2). The data in Table II is generated from the constellation of Figure 3. the 0. nJ4, 2n/4 and 6n/4 rotations in Table I, and 
the real part of the possible symbol product combination of ' $2 Each factor of each partial response in the auto 

35 correlation is constructed by a one of the symbol values times a channel estimate. 

In accordance with the invention, components of the cross correlation terms and auto correlation terms of a group 
0? branch metrics at a symbol instant are calculated once and stored in respective registers 564-586. Selective ones of 
the components of the cross correlation terms and auto correlation terms are combined by multiplexers 590-600 and 
selectors 626-636 to form the cross correlation terms and auto correlation terms. State table control 588 determines 

40 which components from registers 564-586 are combined and whether the components are positive or negative. Outputs 
650, 652 and 654 are the cross correlation terms in equation (2) and outputs 656, 658 and 660 are the autocorrelation 
terms in equation (2). Adder 662 combines the cross correlation terms and auto correlation terms to produce a branch 
metric as output 664, The metric counter is incremented and another branch metric is generated from the components 
stored in registers 564-586. This process continues until all of the branch metrics are generated that are based on the 

45 components stored in registers 564-586. Typically, the branch metrics for all branches into states at a symbol instant 
are calculated based on the components stored in registers 564-586. The channel impulse response in blocks 500, 502 
and 504 are updated as is y n representing the next symbol instant. Updated components are calculated and stored in 
registers 564-586, such as by overwriting. The branch metrics for the updated symbol instant are then computed in the 
above manner. The process is repeated for all symbol instants represented by the received data. Other aspects of the 

50 viterbi decoding, including traceback operations and determining final decoded symbols may be done in any known 
manner. 

Once the calculations are updated and stored in registers 564-586, 64 branch metrics are produced as various 
combinations of the values stored in registers 564-586. By stepping through the 64 controls of state table control 588, 
64 branch metrics are produced at output 664, one at a tim . Calculations are updated for the next symbol instant, and 
55 64 branch metrics at the new symbol instant are calculated. 

Ther are six terms in equation (2). The number of terms in equation (2) is independent of the value of r, the 
number of samples per symbol. Thus, the number of registers 564-586 will remain the same regardless of the number 
of samples per symbol. 

For GSM applications, the circuit of Figur 5 can be simplified. Th simplification is possible because the expected 
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values of received signals take on only two real values. One of the ways the circuit of Figure 5 is simplified is that mul- 
tiplexers 590-600 are not required. 

Table III is a tabulations of the number of operations required to calculate branch metrics on a per symbol basis 
using the Euclidean distance method. L represents the channel length, that is the number of taps in the channel impulse 
5 response, and r represents the number of samples per symbol. 



TABLE III 



w 



15 



Operation 


Adds/Symbol 


Multipies/Symbo! 


Combine complex channel impulse responses 


4t-4 l 


Not applicable 


Calculate error vector 


2*r«4 L 


Not applicable 


Square error vectors 


Not applicable 


2(4) L 


Combine squared error vectors 




Not applicable 


Total 


(6t-H)4 L 


2(4) L 



20 Table IV is a tabulation of the number of operations required to calculate branch metrics on a per symbol 
accordance with the present invention. 



in 



25 



30 



35 



40 



45 



50 



55 




Operation 

Calculate cross 
correlation terms 




Calculate auto 
correlation terms 



Combine rotated auto and 
cross correlation 

terms for 4 L metrics 



L(r-l) 



(r-I) 



L 2 +L 



-L 



TOTAL 



-I 



(M) 



L 2 +L 



4-r-L 



4 r 



-L 



Not applicable 



4t* 



L 2 +L 



TABLE IV 



The total number of operations to complete a metric calculation at one symbol instant for the 4-ary case using each 
method is tabulated in the bottom row of each table. 

As compared to the Euclidean distance method, the computation and power savings of the present invention are 
beneficial up to about eight channel impulse response tape at one sampl per symbol. The savings increase for multiple 
samples per symbol, and thus the present invention becomes more advantageous even with a greater number of taps 
in the channel impulse response. Even where the number of operations are about the same for the two methods com- 
pared in Tables HI and IV, the power savings are greater for the present invention. Multiply operations consume approx- 
imately five times as much power as an odd op ration using the same operands. The power savings are accomplished 
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by performing relative few numbers of multiply operations and performing relative more add operations. 



TABLE V 



Operation 


Euclidean Method 


Reduced Metric Computation Method 




Two Samples per 
Symbol 


One Sample per 
Symbol 


Two Samples per 
Symbol 


One Sample per Sym- 
bol 


Acids 


832 


448 


326 


320 


Multiplies 


128 


128 


48 


24 


Total operations 


960 


576 


374 


344 


Energy Units Con- 
sumed 


1472 


1088 


566 


440 



Table V illustrates the savings in operations and power for the illustrated embodiment. An energy unit is that amount 
of energy consumed in an add operation. The branch metric operation is the most redundant operation in a receiver, 
thus power saved in each operation, here about 50% reduction, is a significant savings. 

The invention is particularly useful in communications systems and equipment employing integrated circuits includ- 
ing this reduced computation technique. Such communication systems and equipment have the advantage of reduced 
complexity of circuitry to accomplish the signal processing. Furthermore, the reduced computation time reduces power 
consumption. 

While the control of multiplexers 590-600 and selectors 626-636 has been described as being controlled by state 
table control 588. the invention is not limited thereto. The control bits can be generated in other ways, such as but not 
limited to a counter and logic circuitry; a control processing unit (CPU), a computer program, or stored in read-only- 
memory (ROM). 

Although the invention has been described as useful in computing branch metrics, the invention is not limited 
thereto. The invention may also be used to calculate the difference between branch or path metrics entering each state. 
The branch or path metric differences cab be used to calculate "soft metrics" for subsequent error correction decoding. 

White the illustrative embodiment has been described as applicable to branch metric calculations for a trellis having 
a constraint length of C, the invention is also applicable to a "reduced constraint" trellis vtterbi decoder. 

While the illustrative embodiment of the invention has been described as useful in digital cellular telephony, the 
invention is not limited thereto. The invention is useful in any digital transmission and detection application employing 
forward error correction such a viterbi decoder. By way of example and not limitation, other applications include reading 
digital data from a disk or CD ROM, and communications between earth stations and satellites. 

While the illustrative embodiment of the invention has not been described as incorporating pipelining, one skilled in 
the art would recognize the enhanced computation efficiency available by utilizing pipelining in the design. Pipelining is 
achieved by initiating computation with a new dataset before completing computation with a previous set of data. The 
more latches used in pipelining, the greater the depth of pipelining. Pipelining causes an initial latency in computation 
time required to fill the pipeline, but maximizes usage of resources. 

Claims 

1. A method for determining a plurality of branch metrics at a symbol instant that are comprised of autocorrelation 
terms and cross correlation terms, including the steps of generating components of the autocorrelation terms, stor- 
ing the components of the autocorrelation terms, generating components of the cross correlation terms, storing the 
components of the cross correlation terms, 

characterized by: 

selecting predetermined ones of the autocorrelation components and the cross correlation components, and 
combining the selected autocorrelation components and cross correlation components in a predetermined 
manner to produce a branch metric. 

2. The method as recited in claim 1 , further characterized the step of: 

repeating the selecting and combining steps to produce another branch m trie. 
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3. An integrated circuit including a branch metric calculation circuit including a circuit (526. 538, 552) for generating 
components of autocorrelation terms, at least one register (576, 578, 580. 582, 584, 586) for receiving and storing 
components of the autocorrelation terms, a circuit (530, 544, 558) for generating components of cross correlation 
terms, at least one register (564, 566, 568, 570, 572, 574) for receiving and storing components of the cross cor- 

s relation terms, 

characterized by: 

a selection circuit (590, 592, 594, 596, 598, 600) for selecting predetermined ones of the autocorrelation com- 
ponents and the cross correlation components stored in said at least one registers (576, 578. 580, 582, 584, 
w 586); 

a selective inverter circuit (626, 628, 630, 632, 634, 636) for selectively inverting predetermined ones of the 
selected auto correlation components and cross correlation components; 

a combining circuit (662) for receiving the selected auto correlation components and cross correlation compo- 
nents and for combining the selected autocorrelation components and cross correlation components to pro- 
fs duce a branch metric. 

4. The integrated circuit as recited in claim 3, further characterized in that the selection circuit includes a look-up table 
(588). 

20 5. The integrated circuit as recited in claim 3 or 4, further characterized in that the combining circuit is an adder (662). 

6. The integrated circuit as recited in any of claims 3-5 further characterized in that the selective inverter circuit (626, 
628. 630, 632, 634, 636) provides as an output a positive or a negative version of the input depending upon a con- 
trol select input (632, 634, 636, 638, 640. 642). 

25 

7. The integrated circuit as recited in any of claims 3-6, further characterized in that the selection circuit selects all of 
the autocorrelation components and all of the cross correlation components. 

8. The integrated circuit as recited in any of claims 3-7, further characterized by at least one rotate circuit (506. 508, 
30 510). 

9. An integrated circuit for having a branch metric calculation circuit including a circuit for generating components of 
autocorrelation terms, a circuit for generating components of cross correlation terms, and a circuit for combining 
autocorrelation components and cross correlation components characterized by: 

35 

at least two registers (576. 578, 580, 582, 584, 586) for receiving and storing components of the autocorrela- 
tion terms; 

at least two registers (564, 566, 568. 570, 572, 574) for receiving and storing components of the cross corre- 
lation terms; 

40 a selection circuit (590, 592, 594, 596, 598, 600) for selecting predetermined ones of the autocorrelation com- 

ponents and the cross correlation components stored in said at (east two registers; 
a selective circuit (626, 628, 630, 632, 634, 636) for selectively inverting predetermined ones of the selected 
auto correlation components and cross correlation components: 

a combining circuit (662) for receiving the selected auto correlation components and cross correlation compo- 
45 nents and for combining the selected autocorrelation components and cross correlation components to pro- 

duce a branch metric. 

1 0. The integrated circuit as recited in claim 9, further characterized in that the selection circuit includes a look-up table 
(588). 

50 

11. The integrated circuit as recited in claim 9 or 10, further characterized in that the combining circuit is an adder 
(662). 

12. The integrated circuit as recited in any of claims 9- 1 1 further characterized in that the selective invert r circuit (626, 
55 628, 630, 632, 634, 636) provides as an output a positive or a negative version of the input depencfing upon a con- 
trol select input (632. 634. 636, 638, 640, 642). 

13. The integrated circuit as recited in any of claim 9-12, further characterized in that th selection circuit selects ail of 
the autocorrelation components and all of the cross correlation components. 
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14. The integrated circuit as recited in any of claims 9-13, further characterized by at least one rotate circuit (506, 508, 
510). 

15. A receiver (14), the receiver (14) including a decoder having a branch metric calculation circuit, including a circuit 
5 (526. 538, 552) for generating components of autocorrelation terms, at least one register (576, 578, 580, 582, 584, 

586) for receiving and storing components of the autocorrelation terms, a circuit (530, 544, 558) for generating 
components of cross correlation terms, at feast one register (564, 566, 568. 570, 574, 586) for receiving and storing 
components of the cross correlation terms, a selection circuit (590. 592, 594, 596, 598. 600) for selecting predeter- 
mined ones of the autocorrelations components and the cross correlation components stored in said at least one 
io registers (576, 578, 580, 582, 584, 586). 
characterized by: 

a selective inverter (626, 628, 630> 632, 634, 636) circuit for selectively inverting predetermined ones of the 
selected auto correlation components and cross correlation components, and 
is a combining circuit (662) for receiving the selected auto correlation components, and cross correlation compo- 

nents and for combining the selected autocorrelation components and cross correlation components to pro- 
duce a branch metric. 

16. A receiver as recited in claim 1 5, further characterized in that the selection circuit includes a look-up table (588). 

20 

17. A receiver as recited in claim 15 or 16, further characterized in that the combining circuit is an adder (662). 

1 8. A receiver as recited in any of claims 15-17 further characterized in that the selective inverter circuit (626, 628, 630, 
632, 634, 636) provides as an output a positive or a negative version of the input depending upon a control select 

25 input (632, 634, 636, 638, 640, 642). 

19. A receiver as recited in any of claims 15-18, further characterized in that the selection circuit selects all of the auto- 
correlation components and all of the cross correlation components. 

30 20. A receiver as recited in any of claims 15-19, further characterized by at least one circuit (506, 508, 51 0). 

21, A receiver (14) as recited in any of claims 15-20, further characterized by a transmitter (12). 



35 
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